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Abstract 

A recently proposed theory for diffusion-limited aggregation (DLA), which mod- 
els this system as a random branched growth process, is reviewed. Like DLA, this 
process is stochastic, and ensemble averaging is needed in order to define multifractal 
dimensions. In an earlier work [T.C. Halsey and M. Leibig, Phys. Rev. A 46, 7793 
(1992)], annealed average dimensions were computed for this model. In this paper, 
we compute the quenched average dimensions, which are expected to apply to typical 
members of the ensemble. We develop a perturbative expansion for the average of the 
logarithm of the multifractal partition function; the leading and sub-leading divergent 
terms in this expansion are then resummed to all orders. The result is that in the 
limit where the number of particles ti ^ oo, the quenched and annealed dimensions 
are identical] however, the attainment of this limit requires enormous values of n. At 
smaller, more realistic values of n, the apparent quenched dimensions differ from the 
annealed dimensions. We interpret these results to mean that while multifractality as 
an ensemble property of random branched growth (and hence of DLA) is quite robust, 
it subtly fails for typical members of the ensemble. 
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I. Introduction 

Many natural growth processes generate branched structures. Probably the most 
celebrated such process is difFusion-limited aggregation (DLA). A simple algorithm for 
this type of growth process, introduced by Witten and Sander, has allowed over ten 
years of numerical exploration of its mysterious properties. Theoretically, there 
has been less progress. A number of studies have, with greater or lesser success, 
used real-space renormalization or variants thereof to study the self-similarity of DLA 
clusters. An alternative approach has been to try to deepen our understanding of 
the scaling structure of DLA clusters by emphasizing the "multifractal" nature of the 
growth of such clusters. 

Recently, a quite broad framework has been proposed to study branched growth 
processes.^-*^-*^) This framework relies upon the dynamics of competition between neigh- 
boring branches in these structures. For DLA, one of the authors of this work has 
proposed a method of determining the universal dynamics underlying this competi- 
tion, a method which allows a priori computation of, e.g., the overall fractal dimension 
of DLA clusters. Results of this approach are in good agreement with numerical 
results for DLA. 

In this work, we turn to the implications of this approach for the multifractal 
properties of DLA clusters. Our principal result is that fluctuations in the ensemble of 
possible DLA clusters seem to play a more important role than previously expected in 
determining these multifractal properties. In fact, while multifractality is well defined 
as an ensemble property of DLA, there appears to be a quite precise sense in which 
typical DLA clusters are not multifractals at all! 

This result should be placed in the context of a number of studies claiming devia- 
tions from perfect fractal or multifractal behavior for DLA. A number of groups have 
reported anomalies in the scaling behavior of the very weak growth sites in DLA.*^^^) 
More recently, Mandelbrot and collaborators have proposed that the radius of gyration 
scaling of large DLA clusters may be more complicated than previously believed. ^^^^ 
However, in our case, the scaling anomalies correspond to the most active, strongly 
growing regions of the cluster, and the form of the anomalies can be computed. 
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Diffusion-limited aggregation 

The Witten- Sander algorithm is based on modeUing growth as an aggregation of 
random walkers. Consider a cluster composed of n particles. To generate the n + 1 
particle cluster, introduce a random walker at a great distance from the cluster. This 
random walker will either escape to infinity without encountering the cluster, or else 
will eventually encounter the cluster. In the Witten-Sander algorithm, the particle 
sticks at the point of its first contact with the cluster, thereby forming the n + 1 
particle cluster. The procedure is then iterated to create arbitrarily large clusters (the 
current record in two dimensions, with off-lattice random walks, is in the neighborhood 
of 10^ particles!) 

Structures generated by this algorithm are highly branched and ramified (Figure 
1). Fractal dimensions, defined by the scaling of the average cluster radius of gyration 
Tg with particle number n, n cc rg , are Df 1.71 in spatial dimensionality d = 2 
and Df Ki 2.49 in spatial dimensionality d = 3. In higher dimensionalities, fractal 
dimensions are given approximately by the phenomenological formula Df ^ (d^ + 
l)/(rf+l).(i6) 

DLA is observed to be a good model for a variety of natural growth processes, 
including electrodcposition, viscous fingering, and solidification. Variants of the 
model have been proposed to account for colloidal aggregation, dielectric breakdown, 
and the fracture of brittle media. 

Note that DLA is a stochastic process, due to the underlying stochastic nature of 
the random walks. Properly speaking, for fixed particle number n there exists a large 
ensemble of clusters, each with a certain probability of appearing. If we (naively) 
suppose that each arriving particle can attach to any of the pre-existing particles, 
then there are (n — 1)! different "genealogically distinct" (i.e., with respect to particle 
attachments) members of this ensemble. 

This implies that quantities such as Df must be defined as ensemble averages, 
since not all possible members of the ensemble will have the same Df. For instance, a 
possible, though unlikely, ensemble member consists of a long straight chain of parti- 
cles, which has Df = 1. Unfortunately, researchers have rarely paid much attention to 
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this aspect of the problem, and few numerical studies even carefully report the method 
of ensemble averaging used. 

Multifractality 

In a particular cluster of n particles, there is a well-defined probability that the 
n + I'st particle will stick to any particular pre-existing particle. Indexing the particles 
by z, we write the probabilities of attachment {pi}', Pi will equal zero for any particle 
inacessible to an approaching random walker. These quantities are very broadly dis- 
tributed, being relatively large at the tips of clusters, where arriving particles are quite 
likely to attach, and extremely small deep inside the "fjords" of a cluster, where ran- 
dom walkers will almost never successfully arrive without having previously contacted 
the cluster. If the distribution of {pi} is multifractal, then we expect that 



where the exponent function a{q), which is defined by this relation, should be asymp- 



is defined in terms of the scaling of moments of {pi} with respect to the length scale 
Tg, and not particle number, defining an exponent function r(g). We use a{q) because 
for our purposes particle number is a more convenient basis for defining the exponent 
function. We expect T(g) = Dfa{q).) 

These exponent functions have well-known interpretations in terms of the dis- 
tribution of {pi}- We will temporarily follow the literature and use T{q) rather than 
cr{q)- Suppose that for any particular a corresponding ai is defined hy pi = {a/rg)°^\ 
where o is the particle diameter, and Vg is the cluster radius of gyration. Then the 
total number of particles with pi in a range da, N{a)da, defines a quantity /(a) via 
N{a) = (a/rg)"-^*^"); /(a) is independent of rg for multifractals.'-^^^ The quantity f{a) 
can be obtained by Legendre transformation of T{q), 



n 




(1.1) 



totically independent of particle number nM^^'^'^^'' (Note that usually multifractality 




(1.2) 



f{a{q)) =qa{q) - r(g). 
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The values of / and a at a particular value of q satisfy the tangent condition 



By definition, /(a) is an intrinsically positive quantity, since N{a) > 1. 

Now we must modify these standard relations to account for the stochastic nature 
of DLA. One approach is to define the exponent function (j(g) by a simple average of the 
moments; by analogy with statistical mechanics, we term this "annealed" averaging, 
and the resulting dimensions a^iq) "annealed" dimensions, 



where the brackets (• • •) denote, here and elsewhere, averaging over the ensemble of 
DLA clusters. We can similarly define T^{q) in terms of the scaling with respect to 
the cluster radius-^ Vg by {Yl^=iPi) = ('^/'^g)"^'^^^^- 

One advantage of this procedure is that the corresponding function f{a) retains 
a simple interpretation. Suppose that the expectation value of the number of particles 
with probabilites pi corresponding to the exponent a is {N{a)). Direct application of 
Legendre transformation to Eq. (1.4) gives {N{a)) oc {a/rg)~f^°'\ where /(a) obeys 
Eqs. (1.2-3), and can be thereby generated from T^(q') = Dfaj\^{q). One difference 
with the non-stochastic case is that, since the expectation value of a number can 
take any non-negative value, f{a) can now be negative. Figure 2 shows characteristic 
f{a) functions for toy models, in the non-stochastic case and for stochastic models 
with annealed averaging. Typically in toy models, negative values of / appear for a 
sufficiently far from the value of a, am, which corresponds to the maximum value of 
/, (am represents the most frequently occuring value of the growth probability.) In 
terms of T^{q) or a"^(g), negative values of / thus correspond to large absolute values 
of q. 

^ Unfortunately, Vg is also a stochastic quantity. The ambiguity this introduces is 
one argument for using cr{q), defined in terms of n, for which there is no such ambiguity. 



n 




(1.4) 
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The great drawback of annealed averaging arises if the quantity il{q) = J^^pf is 
subject to large fluctuations between different ensemble members. Suppose that these 
fluctuations are log-normal, so that the probability V{fl,)dQ that a particular ensemble 
member has J2iPi ^^e range [Cl, f2 + dfl] is 

V{n)dn = Poexp l^ ^^^g^ - logf])^ ^ dlogn. (1.5) 

Clearly, the most likely value of logfi, logfi, is not equal to log{0,)S^'^^ This raises the 
possibility that observation of the scaling function a{q) for a typical DLA cluster may 
lead to a different result than a^iq)- However, Eq. (1.5) shows that if the fluctuations 
of O are log- normal, then (logO) = logQ. This suggests that dimensions crg(g) defined 
through "quenched" averaging, 

n 

(log ^ ) = -CTQ (?) log n, (1.6) 

i=l 

may be closer to the typical result of an observation on, for instance, a numerically 
generated cluster. Of course, it is not necessarily the case that the fiuctuations of for 
DLA are log-normal. However, whatever the nature of these fluctuations, quenched 
averaging reduces the impact of rare ensemble members, and thus should yield results 
closer to the behavior of typical ensemble members. 

If crQ{q) gives results characteristic of a typical member of the ensemble, then the 
corresponding /(ct) cannot include negative values of /, which can only be interpreted 
as ensemble average quantities. ^^^-^ Thus if the fluctuations of the ensemble are strong 
enough to lead to negative values of / in cr^{q), then we also expect that (J^{q) ^ 

This does not imply that / < are necessarily unobservable. Suppose at some 
value of q, f{a{q)) (annealed) is < 0. Then in order to see this value of a, one must 
average over an ensemble of at least N ~ (r^/a)!-^' clusters. This will become difficult 
as sizes increase, but for small absolute values of /, negative values of / should be 
observable. 
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It is convenient to introduce the "partition sum" Z{q,a,n) which is defined for a 
specific cluster by 

n 

The definitions above are equivalent to requiring that the annealed dimensions 
satisfy 

{Z{q,aA{qy,n)) = l, (1.8) 
and that the quenched dimensions satisfy 

{\ogZ{q,aQ{q);n)) = 0, (1.9) 
in the limit of large n. By definition, {Z{q, a; n)) = n'^Z{q, n), so that 

{logZ{q,a;n)) = alogn + F{q,n), (1-10) 

where F is independent of a. 

Branched growth model 

The branched growth model introduced in Ref. 11 is based upon an analysis of the 
dynamics of competition between neighboring branches. In standard off'-lattice DLA 
algorithms, any particle that aggregates to a cluster has a unique "parent" particle, 
to which it attaches. Since no particle has two parents, no loops can develop in 
the structure. The position of any particular particle can be uniquely identified by 
specifying which branch it lies upon, on a succession of increasing length scales. 

Consider a branch point, which is any particle that is parent to more than one 
succeeding particle. At a particular point in the development of the cluster, the two 
(higher numbers are possible, but are both unlikely and irrelevant) sub-branches com- 
ing off of this particle have respectively total masses ni, 77-2 and total growth proba- 
bilities piiP2- These growth probabilities represent the total probability that the next 
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particle will stick to any constituent of that particular sub-branch. We define relative 
masses and relative growth probabilities by 



Pi _ Pi 

X =- 



y 



PI+P2 Pb 



ni + 77,2 rib 

where we have also defined total branch mass and growth probability nj), p^, each made 
up of a contribution from the two sub-branches. 

Simple kinematics shows that, neglecting fluctuations in the numbers of particles 
arriving at this branch, we have*^^^-* 

=x-y. (1.12) 



d log Ub 

If dx/dlogUh can be expressed as a function of x and y alone, then there is a closed 
dynamics for branch competition as a function of these two variables. In Ref. 12, it 
was shown that by averaging over the intrinsic stochasticity of the DLA problem, a 
function g{x,y) can be determined (at least in cZ = 2) such that 

=gix,y), (1.13) 



dlog rib 

with g{x,y) a calculable function of x and y alone. By symmetry, g{x,y) = —g{l — 
x,l — y), so that there is a fixed point of the dynamics at {x, y) = (1/2, 1/2). For the 
computed function g{x, y), this is a hyperbolically unstable fixed point, with one stable 
and one unstable direction. (Actually, this qualitative feature can also be predicted 
on general grounds, without computation). The unstable manifold of the fixed point 
leads to two other (stable) fixed points, at {x,y) = (0,0) and (x, y) = (1,1). Thus, 
sub-branch pairs, which are born in a tip-splitting process, compete in an unstable 
fashion, so that asymptotically one of the two sub-branches possesses virtually all of 
the mass and all of the growth probability of the pair. 

Although the dynamics of a branch pair in the x — y plane is thus deterministic, 
the stochasticity of the DLA problem is retained since the original position of a branch 
pair in this plane, i.e., when ~ 1, and a tip-splitting event creates the pair, is a 
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random function. In fact, sub-branch pairs born in most sections of the x — y plane 
will be quickly driven to one or the other of the stable fixed points, at which one of 
the two sub-branches completely dominates the other. In Eqs. (1.12-13), we see that 
for general values of x, y, the fundamental scale over which x and y will change their 
values is logn^ ~ 1, i.e., while the branch as a whole is still of microscopic size. 

However, if the sub-branch pair is originally created quite close to the central, 
unstable fixed point, then it will remain in the vicinity of that fixed point up to larger 
values of logn;,. Such pairs correspond to the large, relatively equal branch pairs 
seen in a DLA cluster. Thus, the large scale cluster structure is sensitive only to the 
distribution of birth probability in the immediate vicinity of the unstable fixed point. 
This probability distribution will be determined by the microscopic dynamics of tip- 
splitting, which do not recognize any special role of the unstable fixed point. Thus we 
expect this probability distribution to be only slowly varying near that point; it can 
be approximated for our purposes by a constant. 

This criterion, and the specific form of g{x,y), defines the random branched 
growth model. Computations with this model are much simplified by the fact that the 
branch competition trajectories born quite close to the stable manifold, which corre- 
spond to the large-scale structure, will be quickly drawn onto the unstable manifold. 
This manifold can thus be parameterized by x and y as functions of logn^, or of rib- 
Of course, the manifold is symmetric about the unstable fixed point. It is convenient 
to introduce a parameter e, which parameterizes the original distance (when ~ 1) 
of the branch pair from the stable manifold in the x — y plane. If this distance is e'^, 
where is the eigenvalue along the unstable manifold at the unstable fixed point of 
the dynamics defined by Eqs. (1.12-13), it can be shown that x and y are functions of 
the combination eri;,, and that e plays no further role in the dynamics. 

For practical purposes, we must thus specify x{eni,), y{eni,) in order to compute 
with our model. A useful toy model, "Model Z" is displayed in Figure 3. In this model, 
the unstable manifold follows a straight line away from the fixed point until it strikes 
the line a; = (or a; = 1). It then moves vertically into the stable fixed points. The 
angle of the unstable manifold at the fixed point is the one adjustable parameter for 
this type of trajectory. 
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For more realistic purposes, we can use the unstable manifold computed for DLA 
in two dimensions in Ref. 12. This manifold is displayed in Figure 4. 

Summary of results 

In Ref. 11, it was shown that the annealed dimensions (Ta{q) for a branched growth 
model are determined by the simple implicit formula 



where rj = eni, parameterizes position on the unstable manifold described in the above 
discussion. We have also defined the quantity 'ilj{'r]; q, a) by 

y'^iv) v--y{v)r 

This quantity will appear frequently in our results. For model Z or for the computed 
DLA trajectories, Eq. (1.14) leads for g > to an annealed multifractal spectrum with 
the expected properties, such as negative values of f{a) for large values of q. For g < 0, 
annealed multifractal dimensions are typically not well-defined, due to a divergence of 
the left-hand side in Eq. (1.14). 

In this work, we study the structure of (log Z) for the random branched growth 
model, in order to determine the quenched multifractal dimensions for these models. 
We develop a diagrammatic perturbation expansion for this quantity, and we show 
how sets of terms in this expansion can be resummed to all orders. The form of the 
resulting expression is 



(log Z(g, cT, n)) = (^J drjf]" V(^; 9, f^)^ro(g, cr) logn 

+ri(g,a)(l-n-^(«)) + .-.. 



where the exponent A(g) governing the finite-size corrections is positive and calculable. 
Because (log 2^) = crlogn + F(g, n), the power law terms should be independent of 
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cr; we retain an ersatz dependence on a in Eq. (1.14) because this independence is 
difficult to see directly within our perturbative method. In the limit n ^ oo, the 
power-law and constant terms can be disregarded; thus it appears that the criterion 
/ dr] rf~^il){'q\ 5, cr) = is sufficient to determine crg(q'). Since this is the same criterion 
as that for the annealed dimensions (7a{i)i this implies that CTQ^q) = crji,{q). This result 
is somewhat surprising, and also undermines the interpretation of (JQ{q) as a property 
of a typical cluster, since we expect 0-^(9) to include regions of negative /(a). 

This paradox is resolved by a further feature of this result. It concerns the order of 
the limits g — > cx), n ^ 00. Negative values of / appear for large values of g; however, 
limq^oo ^{q) = 0- Thus, for large g, the power-law terms in Eq. (1.16) do contribute 
logarithms to the result, provided that n is not too large. Thus for moderate values 
of n, a different result is obtained for crQ(g), which proves not to include negative 
values of /. At large n, one will eventually cross over to the asympotic regime in 
which (JQ^q) = cr^(g), but this value of n, which we term Uc, grows very quickly with 
g, as riciq) ~ exp(ae^'^), with a and b positive constants, as shown below explicitly in 
section 4 for a specific model, model Z. Thus for practically accessible values of n, we 
expect never to see negative values of /. 

Of course, this resolution also implies that any typical cluster of finite size n 
should probably not be viewed as a true multifractal, as its exponent function crQ{q) 
actually represents different scaling behaviors above and below the value of q for which 
n = nc{q). One conclusion of our study is thus that multifractality should be properly 
viewed as an ensemble property of DLA, and not a property of individual DLA clusters. 

This paper is comprised of five sections and five appendices. In section 2, we 
develop the perturbative expansion of (log 2). In section 3, we show how families of 
terms within this expansion can be resummed to all orders. In section 4, wc discuss the 
implications of these results for the multifractality of DLA clusters, and we compare 
our analytical results with numerical results for random branched growth. In section 
5, we conclude and summarize. In appendix A, we show how to perform sums that 
play the role of "propagators" in our perturbative expansion. In appendix B, we list 
perturbation theory terms, with their diagrammatic expressions, up to fourth order. 
In appendix C, we list some mathematical identities that are useful in the resummation 
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of section 3. In appendix D we show a simpler method of obtaining some of the terms 
in the final expression for (log Z) . Finally, in appendix E, we show how the expansion 
of section 2 can be recast as a perturbation theory in q — 1. 

2. Perturbative expansion for (log.E) 
Formally, we can write 

(logZ(^,a,n))= J]^^^((Z-in, (2.1) 

iV=l 

so that our problem is one of computing expectation values of powers of Z — 1. We 
know that 

n n q 

Z(g,a,n)=n-^pHE(J^, (2.2) 

where we recall that the index i labels the particles. 

In our model, the only sites which are allowed to grow are at the end-points of the 
branching process. Thus we use Eq. (2.2), substituting a sum over these "elementary" 
sub- branches for the sum over particles i in (2.2). Consider a particular elementary 
sub- branch of the cluster, of total growth probability pe and total mass Ue (Figure 5). 
This sub-branch branches off from its sibling at a particular node, indexed by J; the 
total number of particles in the elementary sub-branch and in its sibling taken together 
is nj. If the elementary sub- branch is the weaker (stronger) of the two sub- branches, 
then it has a proportion x^ejrij) (or 1 — x(ejnj)) of the total growth probability of 
the two sub-branches taken together, and a proportion y{ejnj) (or 1 — y{ejnj)) of 
the total mass of the two sub-branches taken together. Here ej is the random variable 
specifying the original state for the branch point indexed by J when the total number 
of descendants of that branch point was nj = 1. We will sometimes refer to the path 
that leads through the stronger of the two sub- branches as the "main branch", the 
weaker sub-branch corresponds to a "side-branching." 

Since Z can be expressed in such a way that x and y always appear together, it 
is convenient to define quantities f±{en) by 
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f-iv) 



f+iv) 



{l-x{r])y 



(2.3) 



The contribution of our particular sub-branch to Z can be written as 



J 



{ne/nY 



(2.4) 



where jjij = ±, depending on whether the stronger or weaker branch is taken at the 
j'th node. The index j indexes the nodes between the root of the cluster, j = 1, and 
the elementary sub-branch, for which we have taken j = J. Note that ni = n, the 
total number of particles in the cluster. Note also that the index j measures position 
upon a particular path from the root to an elementary sub-branch, so its meaning 
depends upon the exact sequence of {/Xj} chosen. Comparing with Eq. (2.2), we see 
that in Eq. (2.4) we are introducing a slightly different cut-off procedure than that 
used in Eq. (2.2), where TT-e = 1. 

Averaging the quantity appearing on the right-hand side in Eq. (2.4) is difficult be- 
cause Uj is actually a function of all e/. with k < j, since either Uj = j/(ej_inj_i)nj_i, 
if at the j'th node we take the weaker of the two sub- branches originating at that node, 
or else rij = (1 — y(ej_inj_i))nj_i, if we take the stronger of the two sub-branches. 
Thus the individual terms in the product appearing in Eq. (2.4) cannot be averaged 
independently of one another. 

Nevertheless, it is instructive to investigate how a single one of these terms av- 
erages. Consider f_{ejnj). The average of this over the random variable Sj is given 



In the introduction, we stated that the initial distance in the x — y plane from the 
unstable manifold is proportional to e'^. In addition, we stated that the probability 



by 




(2.5) 



where Q^. denotes averaging over e^, and p{ej) is the probability distribution for e^. 



14 



distribution with respect to the distance measure de'^ should be constant close to the 
unstable manifold. Since de'^ oc e'^~^de, we see that 



limp(e)=poe^-\ (2.6) 

e— >0 

where po is a constant. (The reader may consult Ref. 11 for a more detailed discussion 
of this point.) Now suppose that rij » 1. The quantities x{r]), y{rf) — > for » 1, 
(see Figures 3 and 4), and for physical values of cr, we thus expect /-(r/) — > for 
7/ » 1. Thus we have 

OO l-OO 

de,p{e,)f.{e,n,) ^ ^ / dr, rf-'f-iv) = ^((/-)), (2-7) 
where we have defined a normalized expectation value ((• • •)) by 

OO 



ifm^ / dvV^-'fiv)- (2i 







This expectation value does not depend on Uj, nor on anything else except the form 
of the unstable manifold. The corrections to Eq. (2.7) are of higher order in nj^. 
Ignoring these corrections, we write 

(/-(e,%)) = (^)((/-)) (2.9) 

where the remaining true expectation value {nj'^) depends only upon random variables 
e/e with k < j. 

The averaging of /_|_ is performed rather diff'erently, because lim^^oo/+(^) = 1- 
Thus 

/•OO poo 

/ dejp{ej)f+{ejnj) = / dejp{ej) [1 + {f+iejUj) - 1)] 
^0 Jo 

pOO 

=1+ / dejp{ej){f+{ejnj) -1) 

, (2.10) 

-l + ^(((/+-l))) 
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where we have used dep{e) = 1, and have introduced 

g+{rj) = U{rj)-l. (2.11) 

It is now natural to expand (log 2) perturbatively in /_ and g+. Formally, we 
can associate a dummy parameter 6 to each factor of /_ or expand order by order 
in 6, and set 5 = 1 at the end of the computation. 

Zeroth order in 6 

At zeroth order in 5, only one elementary sub-branch in the cluster contributes, 
that in the stronger sub- branch at each node, starting at the root (Figure 6). The 
contribution p%/ {rie/nY from this sub- branch is given by 

J J 

\[f+{^jnj) = n (1 +^7+(6,n,)) (2.12) 

where the index j now denotes nodes along this main branch of the cluster. To zeroth 
order in this is simply p%/{ne/nY = 1, so that Z — 1 = 0(5), and to zeroth order 
in 5, (logZ) = 0. 

First order in 5 

To first order in 5, we can approximate 

(logZ)^(Z-l), (2.13) 

because the other terms in the expansion of the logarithm arc at least of 0(^2). The 
terms with only one factor of /_ represent elementary branches removed at some point 
from the main branch by the choice of the weaker branch at only one node. If there are 
J nodes in the main branch, there are J such "first-order" sidebranches. Since we are 
interested in terms of first order in 5 only, for these terms the factors of /+ appearing 
in the partition function can be replaced by 1, because f^ = l + gj^ = l + 0{5). Thus 
the term of 0{5) that is proportional to /_ is 



Pi 
{ue/nY 
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{^ogZ),,^f_=J2(^){{f_)), (2.14) 

where the subscript on the left-hand side represents terms of first order in S that are 
also proportional to ((/-)). The index j represents nodes on the main branch of the 
cluster, starting at the seed (or root) j — l- 

There are also terms of first order in 5 that are proportional to {{g+))- These terms 
include no factors of /_ , and refer to the elementary sub- branch at the end of the main 
branch. The actual partition sum contribution from this sub-branch is represented as 
a product of factors of /+, one at each node. To first order in 5, one chooses one of 
these nodes to contribute a factor of g+, and all others contribute a factor of 1. Thus 

(logZ)i,«,^=^(^\((^+)), (2.15) 

where, again, the index j refers to nodes on the main branch of the cluster. 

Although the term oc g+ is thus identical in form to that oc /_, Eq. (2.14), the 
origin of the two terms is quite different. Equation (2.14) gives the 0{5) partition 
sum contribution of all of the elementary sub-branches along the "main line" of the 
sidebranches of the main branch, while Eq. (2.15) gives the 0{5) contribution from 
the elementary sub-branch at the end of the main branch itself (Figure 7). 

We can also represent these terms graphically. We represent a g.^. "vertex" by an 
open circle, and an /_ vertex by a solid circle. Their sum 

^ = f-+g+ (2.16) 

(from Eq. (1.15)) is indicated by two circles connected by a short vertical line. A set 
of nodes on a branch over which nj'^ (and later more complicated functions of Uj) is 
summed is indicated by a solid horizontal line. The left side of a diagram indicates 
the root, and the right hand side represents structure successively further down the 
branching tree. Thus the (logZ)i terms that we have written in Eqs. (2.14-15) above 
are indicated by the diagrams in Figure 8. 
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The perturbation expansion we are developing is analogous to field theoretic per- 
turbation series, with the factors of /_ and g+ playing the role of vertices, and the 
sums of playing the role of propagators. We shall see below that at higher order, 
both of these objects become more complicated. One difference with field theory is 
that topologically, the diagrams always have the shape of branched trees, with no loops 
in the structure. 

To evaluate these terms, we must evaluate ^7^)' where the sum is along the 
nodes of the main branch. We shall evaluate this sum by a somewhat roundabout 
path. Consider the quantity log(n/ne), where is the number of particles in the 
elementary branch at the end of the main branch. By the definition of the y parameters 
(Eq. (1.10)), we can write the identity 

J 

log(n/ne) = - log J] (1 - yiejuj)) , (2.17) 

where the index j ranges over all of the nodes on the main branch, of which J is 
the last. Taking the expectation value of the right-hand side, and using the methods 
we have developed above for computing these expectation values, we see (referring to 
Eqs. (2.7-8)) that 

logn - logne = - 5^ ( ^ \ ((log(l - 2/))), (2.18) 

so that 

E^^^ + (2.19) 

with A = — ((log(l — y)))~^ and oq = logne/((log(l — y))). The parameter A is thus a 
function only of the unstable manifold in the x — y plane. This is not the case with oq, 
which is cutoff dependent, and thus sensitive to small-scale details of the theory. We 
shall see below that such non-universal constants do not affect results for multifractal 
dimensions. 
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Thus our final result is that 

(logZ)i = ((/_ + (7+))(Alogn + ao). (2.20) 

We should remark that the non-universal constant po, and the exponent no longer 
appear explicitly in our formulae, although they do appear in intermediate steps in this 
computational method, and i/ appears in the definition of the average ((• • •)). These 
quantities do not appear explicitly in final results at any order in S. 

In Eqs. (1.15) and (2.16), we defined the quantity /_+/+ — 1 = f-+g+ = ip, and 
stated that the criterion {{^|;)) = determined the annealed multifractal dimensions 
(^AiQ)- (Recall that like /_ and g^, is a function of q and a.) Thus, in the 
limit ?i oo, our 0{d) result for the quenched multifractal dimensions, determined 
by the criterion (logZ)i = 0, is identical to our result for the annealed multifractal 
dimensions, so that crQ{q) = (J^(g) + 0{5'^). 

Second order in S 

Although similar in spirit to the computation of (log2^)i, the computation of 
(logZ)2 introduces some new complications. Let us start by considering the term 
arising from (Z — l)2^ocg2 • This term arises from the elementary sub-branch at the end 
of the main branch; in this case we are computing the 0{S^) term in its contribution 
to the partition sum. 

Specifically, we must compute 

- l)2,oc4 " iJ2J23+(^''^k)g+{ejnj)). (2.21) 

k>l j>k 

Graphically, this term is represented by a horizontal bar interrupted by two open 
circles, representing the two factors of g^, one downstream from the other on the 
main branch. Note the presence of two horizontal segments, corresponding to the two 
sums in Eq. (2.21) (Figure 9). 

In Eq. (2.21), we start by averaging over e^, thereby obtaining 



19 



{Y.Y.9+i^k''k)9+{ejnj)) = ^^+(efenfe)^)((^+)). (2.22) 

k>lj>k k>lj>k ^ 

Using the result from Eq. (2.19) for ''^J'^ -i we see that 

^) = (Alognfe+i + ao) = (Aflogn^ + log(l - y{eknk))] + ao). (2.23) 

j>k 

where Uk+i refers to the number of particles in the stronger branch at the k'th node 
on the main branch. Then, by summing over j and averaging over e^, we obtain 

(Z-l)2,oc4 -(T.^{^[Miognk + {{g+log{l-y)))]+aoM^ (2.24) 

\ k>l I 

In order to perform the final sum, we must compute logn^/ri^j)- This can 
be done by an iterative procedure based upon our result above for Ylij ''^J'^ -i ^q- (2.19) 
(details are given in appendix A.) The result is 

(E ^^^) = (^2,2 log^ n + A2,i logn + A2,o) , (2.25) 

j>i "^J 

where n = ni is again the total number of particles in the entire cluster, and A2,o is 
a non-universal constant. By contrast, A2,2 and A2,i are functions only of the form of 
the unstable manifold. The generalization of Eq. (2.25) is straightforward 

, N-l N 

(E '° ' = E log" »• (2-26) 

j>i i M=0 

In this notation, all of the {Ajv,o} are non- universal; the first of these is Ai^o = o-q- 
Furthermore, A = Ai^i. In appendix A, an iterative method which can be used to 
compute arbitrary numbers of these coeflBcients (except, of course, for the non-universal 
ones) is demonstrated. For the purposes of this work, we need only A = ((log(l — 
A2,2 = A/2, and A2,i, which is given by 
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((log^(l-y))) _ 
^^'^ - 2((log(l-,)))^ = ^ 

Applying these results to the computation of the right-hand side of Eq. (2.24), we 
obtain 

+ (Alogn + ao)A((g+ log(l - y))){{g+)) . (2.28) 

+ {Xlogn + ao)ao{{g+)f 

Another term, {Z — 1)2, oc/^ > corresponds to partition sum contributions from 
elementary sub-branches with two weak ancestor nodes. In this case, since we wish to 
consider only terms of 0(5'^), we take no factors of g+ whatsoever. The diagrammatic 
representation of this term is shown in Figure 10. A short vertical bar is added before 

the second horizontal line, to indicate that one of the summations of n~'^ takes place off 
of the main branch. The calculation is entirely analogous to that leading to Eq. (2.28), 
and the result is 

{Z - l)2,oc/2 =(A2,2 log' n + A2,i log n + X2,o)X{{f-)f 

+ (Alogn + ao)A((/_ logy)) ((/_)) . (2.29) 
+ (Alogn + ao)ao((/-))' 

The mixed terms (Figure 11), which contain one factor of and one factor of 
/_ (in either order), have the same structure. 



{Z - l)2,oc/_.,+ =(A2,2 log' n + A2,i logn + A2,o)2A((/_))((^+)) 

+ (Alogn + ao)A[((/_ logy)){{g+)) + {{g+log{l - y))){{f-))] • 
+ (Alogn + ao)2ao((/-))((^+)) 

(2.30) 

There are also contributions to (logZ)2 from the —{1/2){{Z — 1)^) term in the 
expansion of the logarithm. These reflect cross-products between two elementary sub- 
branches, with the partition sum contribution of each being taken to 0{5). The di- 
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agrams appearing are shown in Figure 12. Note that symmetry factors appear mul- 
tiplying some of the diagrams. Also, new vertices of the type ((^+)), {{f-)) appear, 
arising from cases in which g+ or /_ factors are taken at the same vertices in the two 
different "replicas" of Z — 1. Such vertices are indicated by adjoining circles, displaced 
perpendicularly. The contribution oc g"^ is 

-^((^-l)')2,oc4 = -(A2,2log'n + A2,ilogn + A2,o)A((^+F 

-(Alogn + ao)A((^+log(l-2/)))((^+)) 

(2.31) 

- (Alogn + ao)ao((t/+))^ 

- ^(Alogn + ao)((^5-))' 

while that oc is 

— {{Z- l)\,ocP_ = -(A2,2 log' n + A2,i logn + A2,o)A((/_F 

- (Alogn + ao)A((/_ log(l - y))){{f-)) 

(2 32) 

- (Alogn + ao)ao((/-))2 
-^(Alogn + ao)((/!)). 

Note that the vertex ((/_ log(l — y))) appears in this term, rather than ((/_ logy)), as 
in Eq. (2.29). Finally, the mixed term is 

— {(Z- l)')2,oc/_..+ = -(A2,2 log' n + A2,i logn + A2,o)2A((/_)) {{g+)) 

- (Alogn + ao)[A((/_ log(l - y))){{g+)) 

+ A((/-))((^+log(l-y)))] (2.33) 

- (Alogn + ao)2ao((/-))((^+)) 

- (Alogn + ao)((/- ■g+)). 

Collecting terms from Eqs. (2.28-33), we finally obtain 
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(logZ) 



A((#((/_log(-^)))-i((V^2) 
± y / 



(Alogn + ao) (2.34) 



where, again, ip — f- + g^. Note that the log n term has cancelled out, as have the 
non-universal contributions to the logn term (oc aoAlogn in the above Eqs. (2.28-33)). 

Structure of 5 expansion 

In appendix B, we continue this expansion, deriving (log 2)3 and (log 2)4. In 
general, we expect that the 0{S^) term in the expansion will have the form 

TV 

(log Z)n=Y1 I^n,m log"^ n, (2.35) 

M=0 

where the coefficients f3N,M will include "universal" (independent of the {Aat.o} of 
Eq. (2.26)) and non-universal terms. It is clear from the structure of the expansion 
that the leading coefficients I3n,n do not contain any non-universal terms, as the 
leading coefficient \n,n in Eq. (2.26) for any sum ^j^og^~^ rij /n^- is universal (sec 
Eq. (A. 8) below). We also note that with the exception of the 0(5) term, we have 
/5jv,Ar = to fourth order. Below we will show that this is true to all orders. The 
situation is summarized in Table I, which shows, as a function of order in 5 and power 
of log n, the terms appearing in this series. 

Unfortunately, this series does not allow us to directly determine the quenched di- 
mensions crQ{q). Recall that these are determined by the requirement that (log .2^) — > 
in the limit n — > 00 (more precisely, that it be bounded above and below by quantities 
which go to neither ±00, see Ref. 20 for a discussion). Since each individual term 
log^ n diverges as n ^ 00, and since there is no obvious relation between the coeffi- 
cients of different powers of log n in the series for (log Z) , we cannot extract unique 
results for crQ{q) from this series. 

For this reason, we turn in section 3 to the resummation to all orders in 5 of the 
leading and sub-leading divergent terms in this series. We shall see that the annoying 
higher order logarithms can be safely resummed to decaying power laws, allowing 
unambiguous determination of CTQiq) in the limit n — > 00. 
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Resummation of "one-propagator" series 

Before turning to the general re-summation of the series, we would like to remark 
on a simpler resummation, which is analogous to loop expansions in ordinary field 
theory. Consider all terms in the series with only one summation over n~'^ . These 
terms can be of any order in 6, as vertices can be multiplied together at the same 
node, as in the 0(5^) terms above. The associated diagrams are shown in Figure 13. 
This "one-propagator" result can be easily derived: 

(logZ)ip= ^ ^ ^ ((V^^))(Alogn + ao). (2.36) 

Ar=i 

We have used the subscript Ip to indicate that this is the one-propagator term. Taking 
the sum inside the brackets (()) yields 

(logZ)ip = ((log(l + V')))(Alogn + ao). (2.37) 

This result has two appealing properties-in the first place it leads to a result for 
CTQiq) in which no negative values of / appear. It is thus an appealing approximate 
formula for the dimensions of a "typical" cluster. Also, as we shall see in section 4, 
the ctq (q) resulting from setting (log Z)ip = is quite close to numerical results in the 
case of our toy model, model Z. 

However, if one tries to extend this approach by resumming the "two-propagator" 
terms, and so forth, one encounters the same problem as with the (5-series, that higher 
orders in logn are also generated. Thus this expansion technique suffers from the 
same difficulty as the d expansion, that unique values of c"Q(g) are impossible to 
obtain without resummation. 

3. Resummation of the 5-series 

We have shown that the form of the 5-series is: 

oo QO N 

(logZ) = ^''{^os2)n = E E /5w>Mlog^n, (3.1) 

N=l N=l M=0 
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where each (3n,m is of 0{5^). Our procedure is now to resum this series by summing 
first the most divergent terms M = N at all orders in 5, then the second-most divergent 
terms M = N—1, and so on. Introducing yet another subscript, £ (leading), to indicate 
the most divergent sum, and the subscript s£ (sub- leading), to indicate the next most 
divergent sum, we have 

oo 

(logZ),= Y,pN,Nlog''n, (3.2) 

N=l 

and 

oo 

{log Z) si = ^ (3n,n-i log^~^ n. (3.3) 

JV=2 

In principle, we could continue with sub-sub leading terms, but in practice we 
will restrict ourselves to computing only these two terms. Note that our resummation 
procedure corresponds to summing down the diagonals of Table I. 

Summation of leading terms 

The leading logarithm at each order in S will be universal, as pointed out in section 
2 above. Let us start by computing {Z — 1)^. In order to compute this, we will need 
the result from appendix A: Xn,n = A/A^, where the reader might recall the definition 

of Xn,m 

1 N-l N 

(T, " t "^ >=EVMlog"r.. (3.4) 

j>l M=0 

Here the sum is along a horizontal bar of some diagram. Since we are taking the 
leading order in logn in the entire diagram, we wish to take the leading order in each 
propagator. No factors of 1 — y or y will appear in the leading order, since each such 
factor would take the place of a factor of log n. Also, we take the leading logarithm of 
each term of the form of Eq. (3.4), which will have as a coefficient one of the {Ajv.at}- 
Consider a diagram contributing to {Z — 1)^ consisting of n\ factors of /_ and 
n2 factors of g^. There will be a total of A = ni + n2 horizontal bars, or propagator 
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sums, in such a diagram (see Figure 14). If we call the contribution of this diagram C, 
then 

C=j^{{f-rH{9+)r^og''n, (3.5) 

because /N\ = Hiv'^i ^n',n'- 

Since any vertex can be taken to be either /_ or g+ without changing other 
features of the term, we immediately obtain 

oo , jv 

(Z-l),= J]— ((V^))^log^n = exp(A((#logn)-l. (3.6) 

N=l 

Now let us turn to — 1/2((Z — 1)^)^ from the expansion of (logZ) in Eq. (2.1). For 
the sake of argument, let us consider terms in which each factor of Z — 1 has at least 
one factor of /_ . We will say that the first factor of Z — 1 has h factors of g+ before 
the first appearance of /_, and nii factors of either /_ or after the first appearance 
of /_. The second factor of Z — 1 has I2 factors of g^ before the first factor of /_, and 
777,2 factors of either /_ or g+ after the first factor of /_ (see Figure 15). 

Now we average these two factors together. Note that two factors of /_ appearing 
in the main line cannot appear at the same node; were they to be averaged together, 
the diagram overall would lose one propagator, and thus one power of log n, without 
losing a factor of 5, which counts the total number of factors of /_ and y^. Such 
diagrams first contribute at sub-leading order. Similarly, two factors of g+ cannot 
appear at the same node. 

Suppose that the factor of /_ coming from the second Z — 1 is further down 
the main line than that coming from the first Z — 1. Then in any diagram resulting 
from this product, there will be a number Z < ^2 of g^ vertices between the two /_ 
vertices, and L — l to the left (further up the main branch) of the first /_ vertex, where 
L = I1 + I2 (see Figure 16). Now we must consider the vertices off the main line. The 
structure of the diagram does not depend on whether these are /_ or g+ vertices, since 
we average no factors of log(j/) or log(l — y) with these vertices at leading order. Thus 
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we will sum together these two types of off-main line vertices, and attach mi vertices 
ip to the first factor of /_ , and m2 such vertices to the second factor of /_ . 

In order to compute the diagram, we must keep track of the number of factors 
of logn appearing for each propagator; at leading order we take only the propagator 
terms multiplying one of elements of the set {Ajv at}. The total number of propagators 
to the right (on the main line) of the "upstream" factor of /_ is m2 + / + 1, while 
the total number downstream of the first factor of g+ to the left of this /_ vertex is 
mi + 1712 + / + 2. (Note that we use "upstream" to mean closer to the root, not closer 
to the elementary sub- branches.) A simple computation, keeping careful track of these 
factors, gives the contribution Wi of this particular diagram as: 



where M = mi + m2 and we defined above L = li + I2. 

Of course, the two original factors of Z — 1 can be multiplied together in various 
distinguishable ways to make this diagram, corresponding to the different origins of 
the L — I factors of g+ on the left hand side of the diagram. Since li of these come 
from the first factor of Z — 1, and L — I — li come from the second factor oi Z — 1, the 
total number of relative permutations Np is 



Finally, to obtain the expectation value of the product of the factor Z — 1 con- 
taining li factors of Qj^ in the main line and mi factors of ip off the main line with the 
other factor oi Z — 1, which contains I2 factors of on the main line and m2 factors 
of '4> off the main line, we must sum over I < I2 and permute the choice of whose main 
line factor of /_ is upstream (see Figure 16). If we are interested in contributions to 
— (1/2)((Z — 1)^) with at least two factors of ((/-)), we must then sum over mi, m2, 
Zi, and l2- This yields 



1 1 (M + Z + 1)! 

(m2 + Z + 1)! mj (M + L + 2)! 

■ {{f-rig+rm^iwogn) 



M+L+2 



(3.7) 




(3.8) 
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OO OO CXD OO ^2 



■i((2-i)W-F=- E E EEE 

mi=0 m2=0 Zi=0 ^2=0 «=0 

'L-l\ 1 1 (M + / + 1)! (3-9) 



^1 J (m2 + / + l)!mi! (M + L + 2)! 
•((/-f((^+))^#r(Alogn)^+^+^ 

In addition, the term including only one factor of /_ is easily written as 



OO 00 00 ^2 



'L-A 1 1 {rm + iy. (3-10) 



Zi J mil l\ (mi + L + 1)! 
((/-»((^+r((V'rnAlogn)-^+^+\ 



while that with no factors of /_ at all is 



1 1 / 7" \ 1 

--{{z - i)^),,oc((/_)>o = - 2 E E / Ty((^+))^(^i«g^)^ (3-11) 

The sums are tedious but elementary. In appendix C we illustrate the performance of 
the types of sums appearing in this work. Performing the sums and adding Eqs. (3.9- 
11) we obtain the quite simple result, 

~{{Z - ir)e = -\ [exp(AW logn) - 1]^ = -\[{Z - (3.12) 

Remarkably, at leading order, the process of taking the expectation value of 2^ — 1 
commutes with taking the square of this quantity. We will now address the origin of 
this identity. 

Let us first consider a simpler case, with only one factor of Z — 1. The leading 
order X> of a diagram with I factors of gj^ on the main line, followed by one factor of 
/_, followed by m factors of i/^, can be written as 
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D = Ar(A((/_))logn)(A((^+))logn)^(A#))logn)-, (3.13) 

where the factor M = l/(m + / + 1)!. The factor N can be written as a nested integral 
over / + 1 variables x^: 

with xq = 1. In this product, the I integrals on the left correspond to gj^ vertices, 
while the last integral corresponds to the /_ factor and the accompanying insertion of 
m additional factors of logn associated with the ij) vertices. 

Now consider a term arising from the product of a factor of 2^ — 1 containing li 
factors of Qj^ in the main line and mi factors of if) off the main line with another factor 
of 2 — 1 containing I2 factors of (7_|_ and m2 factors of if). Again, we choose L = li + I2 
and M = mi + 1712- Each factor of 2 — 1 in addition has one /_ vertex on the main 
line. At leading order, any diagram T> resulting from this product will have the form 

V = X{X{{f-)) \ogn)\X{{g+)) logn)^(A(W) logn)^, (3.15) 



where the combinatorial factor A/" depends on the precise way in which the diagrams 
are multiplied together. It is simple to write a formula for A/". Let us assign variables 
x^^"^ to the main line vertices arising from the first factor of — 1, and variables x^^^ 
to main line vertices arising from the second factor oi Z — 1. Then the factor J\f for 
any particular diagram can be written analogously to Eq. (3.14). For instance, if all 
of the factors of g+ coming from the second factor oi Z — 1 are further down the main 
line (to the right of) the factor of /_ arising from the first factor oi Z — 1, then we 
immediately see that 
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^1+^ mi! 



= n/ 



(M + /2 + 1)! 1 



' rlJ'A ! %/..,+/^'-+^^ (3-16) 

777,2! 



(M + L + 2)!mi! (m2 + /2 + l)!' 

with Xq^"* = and Xq"*^^ = 1. Now the aUowed permutations of the vertices corre- 

spond to all permutations of orderings of {.x^^'*}, {x^^} within the multiple integral 
that preserve the orderings 1 > x^^"* > x')2^ ■ ■ ■ and 1 > x^^'^ > x'^2'' ' ' But this allows 
one to factor the two nested integrals over the points {x^^^} and {x*^^^} completely, so 
that the sum over permutations P oi N can be expressed as 



^1+^ mi! 

1 1 



(mi +/i + 1)! (m2 + /2 + 1)!^ 



where x^^ = 1, which deconvolves the two nested integrals, allowing their factorization 
into terms of the form (3.14). This provides a direct proof of the identity (3.12). 
Using the same argument, we can easily show that at leading order, 

{{Z - = [{Z - 1),]P = [exp(A((V')) logn) - 1]^ (3.18) 
so that the leading order contribution to (log Z) is 

(3.19) 

= A ((V')) logn. 
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We therefore arrive at the result 

(logZ), = log(l + (Z-l),). (3.20) 

We shall see that this is true also at sub-leading order, provided we restrict ourselves 
to terms diverging logarithmically with n. 

Summation of sub-leading terms: Universal 

Sub-leading terms are those with one less power of log n than 5. There are two 
ways in which such terms can arise: 

1. In a term arising from {{Z — with p > 1, more than one factor of /_ or (7_(_ 
may appear at the same vertex. Since these two or more factors are integrated 
together, the net result is that the order of the term in logn is reduced with 
respect to the order in S. If more than two factors of /_ or appear at the 
same vertex, or if more than one vertex contains more than one factor, then at 
least two factors of log n will be lost with respect to factors of 5, so the term will 
be at most of sub-sub-leading order. Thus we need only consider diagrams with 
one vertex possessing two factors of /_ or with all other vertices containing 
at most one factor. 

2. Recall that the form of the propagator sum along a particular sub-branch up to 
a given vertex is ^^=0 ^n,m log^^ n\ where n' is the number of particles below 
the vertex in question, either on the weak or the strong side (See Eq. (2.26)). 
Thus either n' = y{en)n, or n' = (1 — y{en))n, where e is the random variable at 
the vertex being explicitly considered, and n is the total number of descendants 
below that vertex. There are two ways in which the order of the propagator in 
logn may be reduced by one. The first way is for a factor of logy or log(l — y) 
to be taken in the highest order in logn' term (M = N) in the propagator. The 
second way is for the second highest term in log n' (M = A?" — 1) to be taken, but 
without including any factors of log(l — y) or logy. 

Terms of type 1 

An example of a term of type 1. is provided by (logZ)^^ g^02^^. There is no 
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contribution from {Z — 1) to this term; the first contribution arises from — (1/2)((Z — 
1)^) (see Figure 17). It is easy to see that 



As in the above, M = mi + m2, and L = li + I2. The sums are straightforward, and 
yield 



{{Z - l)')siM{fl)) = JFTJ^ [exp(2A W logn) - exp(2A((^+)) logn)] . (3.22) 



The sub-leading terms proportional to (((7^)) and ((/_(7_|_)) are computed in exactly 
the same way, and combine with the term proportional to ((/£)) to yield 



{{Z - lf)se,ocm = [exp(2A((# logn) - exp(2A((^+)) logn)] . (3.23) 



Now we must compute the sub-leading contribution oc ((V'^)) arising from {{Z — 
1)P), with p > 2. A typical diagram is shown in Figure 18. There are a total of 
p{p — l)/2 choices of the two factors oi Z — 1 which contribute the factors of /_, 
fif_(_ which will be integrated together. The computation is considerably simplified by 
the fact that the remaining factors of ^ — 1 can be averaged separately of these two 
factors, by a simple extension of identity (3.18), which we used above to factor the 
leading order terms. Thus 




00 00 00 00 



((^ - i)^>s^,oc((/2» = 

li=0 l2=0 mi=0 m2=( 



(3.21) 




{{Z - If) s£,oc{{^^)) 



pip - 1) 
2 



{{Z - lf)se^^^(^^2))[{Z - l)e] 



(3.24) 



leading to 
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(3.25) 

Using Eqs. (3.6) and (3.23), and substituting from Eq. (3.25), we then obtain 

aog2)sAoc((v>^» = - [1 - exp(-2A((/_)) logn)] 



1 



_ (^1 n J . 



(3.26) 



4((/. 

Recall the definition of ((/-)): 

where x{ri) and |/(?7) refer to the unstable manifold for branch competition. Clearly 
((/_)) > 0, which implies that the correction to (logZ) displayed in Eq. (3.26) ap- 
proaches a constant as n — > oo. Thus this correction will not affect the values of the 
multifractal exponents crQ{q) (defined in Eq. (1.6)), provided that they are computed 
in this limit, as the terms logarithmic in n in (log 2^) will still dominate. 

Terms of type 2 

Now we turn to sub-leading terms of type 2, for which a sub-dominant term in 
one of the propagators is taken. We first consider such terms involving factors of 
log(l — y) or logy (to be averaged with /_ or ^+) arising from the expansion of a 
leading propagator term, 

Aiv.iv {log[(l - y)n]}'' = AAr,Ar[log^ n + iVlog(l - y) log^"^ n + • • •], (3.28) 

with a similar result for logy. Since Xn,n = ^/N, the log(l — y) in Eq. (3.28) multiplies 
Alog^""^ n, which is comparable to a leading propagator term with one less vertex. 
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We will first consider terms containing the factor ((/_ log y)) . Looking at such a 
term arising from {{Z — 1)^), we see that the different factors oi Z — 1 interact only- 
through the various possible orderings of their and one /_ vertices on the main line; 
thus we can factor these diagrams, and consider only the contribution from a single 
factor of Z — 1 (see Figure 19). This is 



1=0 m=l ^ ' )• 

oo oo 

+ E E 7^^(A((/-log|/)))(A((/_))logn)(A((^+))logn)' 
•(A((V'))logn)— ^ 

(3.29) 

The first term on the right-hand side of Eq. (3.29) corresponds to the case where the 
factor of /_ on the main line is averaged with log y, and the second term corresponds 
to the case where a factor off of the main line is averaged with logy. The sums yield 



{2 - l)seM{f- logy)) = A^((V')) log^((/- logy)) exp(A((V')) logn). (3.30) 

Turning to {{Z — 1)^), we see that there are p choices of which factor of Z — 1 
contributes the sub- leading term oc ((/_ log y)) . Thus 



(log2)s£,cx((/_ log J/)) = - ^)s£M{f- logy)) E Z P^(^ " 1)^]^"^ 

p=l ^ 

1 \ (3.31) 



- - l)s^,cx({/_ logy)) (^2+ (Z- 1)^) 

= A2((#((/_ logy)) logn. 



By contrast, terms oc ((/_ log(l — y))) can arise only from a combination of two or 
more factors oi Z — 1, since a factor of /_ is always associated with a weak branch in 
a single partition function, and hence potentially only with a factor of logy. Consider 
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a term arising from a product of p factors of 2^ — 1, each containing U factors of 
in the main hne, followed by one factor of /_ and then rrii factors of ij) off of the 
main line. Let us further suppose that the factor of ((/_ log(l — j/))) arises from the 
factor of /_ appearing in the main line of the i = 1 factor. Then writing L = li 
and M = ^j^rrii, we see that any particular sub-leading contribution proportional to 
((/_ log(l — yf) (which we term T>) will be given by 



V = AfiX{{g+)) lognnxm logn)^(A((/_ log(l - y)))) (A((/_)) logn)^-\ (3.32) 

where the factor jV is determined by the precise ordering of vertices. In fact, it is easy 
to see that will be given by a version of Eq. (3.16) appropriate to p factors of Z — 1, 
with the one addition that there should be a factor of d/dx^^^_^-^ inside the integral over 
the variable (and to the right the factor of (x|^^^j^)™'i/mi!) ). Now the integral of 

all factors to the right of this new operator will yield some number times raised 
to some power; thus the d / dxi^^_^-^ operator can be viewed as acting on only one of the 
factors of 2^ — 1 at a time. It follows that the sum over permutations of J\f can be 
represented as a sum of averages involving two factors of Z — 1 only, with the other 
factors of Z — 1 averaged separately. More precisely. 



(log2)s^,oc((/_log(l-j/))) =X] ^^9 ^ (~^)^ ^(^ - l)^)s£,oc((/_log(l-j/)))[(^ - 1)^]^ ^• 

(3.33) 



2 

p=2 



This naturally represents a considerable simplification of the problem. 

Now we must compute {{Z — l)^)s^,oc«/_ iog(i-y)»- This is given by (see Figure 

20) 



(M-1)!,^,, (3.34) 
.A__^(A((,4)logn)- 

•(A((#logn)^(A((/_log(l-y)))). 
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where the factor of 2 accounts for the distinguishabihty of the contributions of the two 
partition functions. As in the above, M = mi + 777-2, and L = li + I2. The sums may 
be performed directly, and yield 



((/-)) (3.35) 
[exp(2A((V')) logn) - exp(2A((^+)) logn)] . 



Since 



i:,-i)-(p-i)K.-i),r-=-(,^^) ^^^^^ 

= -exp(-2A((V'))logn), 

our result is 



(log2)s^,oc«/_ iog(i-y))> = - A((/_ log(l - y)))^(^j^ [1 - exp(-2A((/_)) logn)] 

= -A((/_log(l-,)))^(l-n-^«/-))), 

(3.37) 

i.e., a constant plus a decaying power law in n. 

Terms proportional to {{g^ log(l — y))) appear at sub-leading order both from one 
factor oi Z — 1 and from two factors averaged together. The first case arises from the 
fact that a g+ factor has a strong branch associated with it downstream, which can 
provide a factor of log(l — y). The combinatorics, and the result, are thus the same as 
for the term oc ((/_ logy)), given in Eqs. (3.29-31). 

The second case involves the averaging of a g+ factor arising from one partition 
function with a factor of log(l — y) arising from the second partition function, leading 
to the same combinatorics, and result, as for the term oc ((/_ log(l — y))) given in 
Eqs. (3.33-37). Thus one obtains 

(log ^),£,oc((g+ log(l-2/))) ^H9+ log(l - y))) {{lp)) 

■ i^^''^'' ~ " exp(-2A((/_)) logn)] 
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We have been discussing sub-leading contributions in which a sub-leading prop- 
agator term is taken within the diagram, leading to factors logy or log(l — y) being 
averaged with /_ or g+. In addition, there are sub- leading contributions in which 
such factors do not appear, but in which Xn,n-i appears rather than Xn,n = X/N. 
In appendix B we show that 

_ l ((log^(l-y))) _ 
2 ((log(l - 

is independent of N. Comparing the expansion of the highest order propagator term, 
Eq. (3.28), with the form of the full propagator, Eq. (3.4), we note that the inde- 
pendence of of Xn^n-i = X' allows us to directly translate our previous results to 
obtain the terms oc A', since the two terms in the iV'th order propagator oc log"^"^ n 
(one including a factor of either logy or log(l — y), and the other a factor of A') 
are similar in form. In particular, we substitute A((/_ logy)) A' ((/_)) in Eq. (3.31), 
A((/_ log(l-y))) ^ A' ((/_)) in Eq. (3.37), and A((^+ log(l-y))) ^ A'((^+)) in Eq. (3.38), 
to obtain 



(logZ)s^,ocA' = X' 



(3.40) 



which is again a combination of a logarithmic divergence, a constant, and a decaying 
power law. 

Summation of sub-leading terms: Non-universal 

At sub-leading order, there are also non-universal terms. These arise from taking 
one factor of ao = Ai,o instead of Alogn in the first propagator sum in a diagram; 
taking non-universal terms in any propagator sum upstream of this will lead at most 
to sub-sub- leading terms. Reviewing the derivation of Eq. (3.6), we see that 

oo ^ 

{Z - l)si,ocao = J^iMi'l')) log^)"^ = «o W exp(A((V')) logn). (3.41) 

7V=0 

Note that this summation starts at N = 0. It follows immediately, using the methods 
developed above to prove factorization of averages of products of Z — 1, that 
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(}ogZ)u,o.a, ^{Z- l),,,ocao Y.^-lY~^[{Z - l),]^"! = ao((#. (3.42) 

Summary of leading and sub-leading orders 

To summarize, we have the following result for (log 2), combining the leading 
order, sub- leading universal, and sub- leading non- universal computations (Eqs. (3.19), 
(3.26), (3.31), (3.37-8), (3.40), and (3.42)): 



2((/-) 
1 „ 



(logZ) =\m logn [1 + A {{{g+ log(l - y))) + ((/_ logy))) + A' 

1 [l_^-2A((/_)> 

2..')) + AW((^log(l-2/))) + A'(W)2 



(3.43) 



where the 0(5"^ logn) terms are sub-sub-leading (see Table I). 

Examining in detail the re-summations performed above, we see that all of the 
terms contributing to the coefficient of logn in Eq. (3.43) arise from separate averages 
of single factors of Z — 1; averages together of two factors contribute to the second 
term on the right-hand side of Eq. (3.43). Anticipating that this is true to all orders, 
we see that^ 

hm ^ = hm i?^. (3.44) 

n^oo logn n^oo logn 

Thus, if we are concerned only with the term proportional to logn, which should 
dominate as n ^ oo, this can be obtained from \og{Z). This quantity is computed by 
much less elaborate means in appendix D. 



^ Using the terminology of statistical mechanics, averages that can be reduced to 
individual averages {Z—1) may be termed disconnected, while those involving averages 
together of more than one factor of Z — 1 may be termed connected. Thus Eq. (3.44) 
corresponds to a sum of disconnected averages only. 
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There are two methods of checking Eq. (3.43), which resums aU terms of the form 
6^ log^ n or log'^"^ n, with N >1. The power laws in n may be converted into an 
expansion in log n, and then compared directly with the S expansion, as developed to 
second order in section 2 above, and to fourth order in appendix B. The reader may 
easily satisfy himself that these two series agree, provided that comparison is confined 
to terms proportional to 5^ log^ n or 5^ log"^"^ n, in this case with 1 < N < 4. 

A second method of checking this result is to recall that the dependence of the 
partition function Z upon a is quite trivial: Z n" (see Eqs. (1.7) and (1.10)). 
However, in our perturbative method this cr-dependence is distributed throughout the 
branched tree, with a factor of y'^ or (1 — yY corresponding to this u dependence 
appearing at each vertex. This considerably obscures the cr dependence, but may also 
be used to check the final result for (logZ), Eq. (3.43). Let us write again Eq. (1.10), 

(logZ(g, (t; n)) — crlogn + F(g, n), (3.45) 

where F has no dependence upon cr. This implies that 

d(logZ) 



da 

and 



= logn, (3.46) 



^^^^ = 0, r>l. (3.47) 
da"- ^ ' 

These relations provide powerful constraints on our result. To check these identities 
against our result for (logZ), Eq. (3.43), it is necessary to use (see Eqs. (2.3) and 
(2.11)) 



da 

and 



^^± = -(^+ + l)log(l-|/), (3.48) 



^ = -/-log(y). (3.49) 
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Note that due to Eq. (3.48), taking the derivative with respect to a effectively reduces 
the order in 8 of terms involving g^^ by one. Thus we can only compute derivatives 
with respect to u of Eq. (3.43) to OiS). 
Let us rewrite Eq. (3.43) as 

(logZ) = A{<i, a) logn + S(?, - n-^^^^-^) + C^q^a) + • • • , (3.50) 

where the coefficients A(g,(j), B{q,a), and C{q,a) can be read from Eq. (3.43), and 
the coefficient C — aodi^)) is the only non-universal one. Taking the derivatives with 
respect to cr, we see after some computation that indeed 

^-^^ = 1 + 0(6% (3.51) 
da 

and 

= 0{5^). (3.52) 

da 

Thus Eqs. (3.46-7) are satisfied provided that the non-universal constant oq = 0. 
It seems that we have succeeded in fixing this quantity. However, we remind the 
reader that many cut-off procedures will effectively break the invariance expressed by 
Eqs. (3.45-47) by introducing a dependence on the number of particles in an elementary 
branch ne, so the result oq = should not be taken too seriously. 

4. Interpretation of results 

In section 3, we re-summed the leading and sub-leading divergent scries to all 
orders in perturbation theory. We found a remarkable result. The resummed series 
consists of three types of terms. In the first place, there are terms diverging loga- 
rithmically with n. These terms are all universal in form. In the second place, there 
are constant terms-at sub-leading order these include non-universal constants. Fi- 
nally, there are power law corrections to (logZ), which involve universal coefficients 
multiplying n~^^ with 
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A = 2A((/_)), 



(4.1) 



a positive quantity. Thus as n oo we expect only the terms oc log n to be relevant 
in determining the quenched dimensions (JQ{q)- Recalling the coefficient of this term 
from Eq. (3.43), we conclude that 



which is equivalent to {{i^iq, o'Q{q)))) = 0. (Setting the other factor equal to zero does 



that determines the annealed multifractal dimensions cr^(q')! This is unsurprising, as 
the logn term in Eq. (3.43) resulted only from products of single averages of {Z — 1), 
as discussed in section 3 or in appendix D. This implies that the true n — > oo quenched 
multifractal dimensions are identical to the annealed dimensions, negative values of / 
and all. 

However, this conclusion is based upon an overly naive interpretation of Eq. (3.43), 
overlooking crossover effects, which we now discuss. Consider again the power law term 
(X with A given by Eq. (4.1). We can understand the role played by this term 

by considering model Z, in which the unstable manifold for branch competition in 
the X — y plane is simply a straight line emanating from the unstable fixed point at 
(1/2, 1/2), which turns vertical upon reaching the a; = or a; = 1 lines (see Figure 3). 
As a function of 77, we can write*^-*^-*^^ 

^ Numerical evaluation of the second factor on the right-hand side of Eq. (4.2) for, 
e.g., model Z (see section 1 or Eqs. (4.3-4)) quickly leads one to the conclusion that 
setting this factor equal to zero does not lead to a physical solution for (TQ{q), or, 
indeed, to any solution for some values of q. 



#(9,^2(9)))) l + X[{{g+{q^aQ{q))\og{l-y))) + ((/_(?, (9)) logy))] 




not lead to an acceptable solution.) 



^ But this is precisely the same as the criterion 
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^=)W-ri% V<1; (4.3) 
'0, ?7 > 1 . ^ ^ 



and 



y U/ry, r/>l. ^^"^^ 

with?7 = [iy/2(l + iy)]. 

Recalling the definition of ((/-)), 

we see that in general, ((/_)) is a function of both q and a. This presents us with 
a paradox, because any exponent A of a power law correction to (log 2^) must be 
independent of cr, due to the invariances expressed by Eqs. (3.45-47). We resolve this 
paradox by realizing that Eq. (4.1) is merely the lowest order in 5 expression for the 
exponent A; we should write (4.1) more precisely as 



A(g) = 2A((/_)) + 0{S') (4.6) 

where A is now the exact exponent of the leading power law correction to (log Z). The 
higher order corrections on the right-hand side of Eq. (4.6) will cancel the dependence 
upon cr arising from ((/_)). 

Let us first compute compute ((/_)) in model Z for cr = 0, 

obtaining 

which goes to zero exponentially in q for large q. 

Now consider (logZ) from Eq. (3.43). For convenience, we set the non-universal 
constant ao = 0. If limg^oo ^{o) = 0, we see that (returning to the case of general cr). 
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lim lim 

g— >oo n— >oo 



(log Z)iij^si 

logn 



lim \m 1 + A(((^+ log(l - y))) + ((/_ log y))) 



(4.9) 



while 



(log 2) 



lim lim 

n— »oo q— »oo log Tl 



= lim X{{^)) 1 + A {{{g+ log(l - y))) + ((/_ log y))) + X' 



-A 



: lim 

q— >oo 



^)) + A((#((V^log(l-y))) + A'((#^ 



A((V^)) +A^ ((#((/_ log(^)))- 



A 



because 



lim lim (1 - n-^(«)) = 2A((/_)) logn + 0(6^). 



n— »oo q— »oo 



(4.10) 
(4.11) 



Thus, for large values of q, the apparent quenched dimensions, which we term a-p{q), 
will be determined by 



A((V^)) + A2 ((#((/_ log(^))) - \{{^'))+0{5^) = 0. 



(4.12) 



These apparent dimensions are the ones that will be seen, e.g., in numerical studies. 
Note that this is simply the naive (non-resummed) result for the dimensions cq^q) 
which cam be obtained from (logZ) = (logZ)i + (logZ)2 + 0(5^) (see Eqs. (2.20) and 
(2.34)). 

Of course, for sufficiently large n, we will always see a crossover back to the true 
quenched dimensions. The crossover value of n for which this occurs, ndq), will be 
determined by 



A{q) lognM = [2A((/_)) + 0(6^)] logn^iq) = 1, 



(4.13) 



or, from Eq. (4.7), 
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ne(?) = exp[^^^;il^ + 0(5°)]. (4.14) 

(Note that the first term in the exponential in Eq. (4.14) is formally 0{5~^).) Thus, 
in practice, the true quenched dimensions will not be observable, even for moderate 
values of q, except for monstrously large systems. Although we have taken an explicit 
form for A appropriate for model Z, we expect the same qualitative results to hold for 
any branched growth model, and in particular, for DLA. 

Suppose that we had not set cr = at the outset of this calculation. It is simple 
to compute ((/-(?, cr))) for model Z, 

«/-(.. <^)» = (i)'"' (l + + 



and 



Um «/_(,,.))) = ^(i)'"' (4.16) 



which is a simple modification of our previous result, Eq. (4.7). For fixed cr, Eq. (4.16) 
leads to the same super-exponential dependence of Uc on q displayed in Eq. (4.14). We 
expect the annealed cr_4(g) to approach a constant as — > oo, while (T-p{q) Pq,P ~ 
l — uKliorq^ oo.*^^^^ It is not, however, legitimate to introduce these functions a{q) 
into Eq. (4.16), because once a becomes large, we must properly include the effects of 
the next order term in 5 in Eq. (4.14), which will offset this spurious u-dependence. 

Numerical results 

In addition to the analytical computation of (log 2), to which the bulk of this 
work is devoted, we have also computed this quantity, for the random branched growth 
model, by a Monte Carlo method. For these purposes, it is convenient to use model 
Z, which can be defined (see Figure 3) by writing x and y as functions of rj as in 
Eqs. (4.3-4). The curve in the x — y plane parameterized by these functions (and their 
complements 1 — x, 1 — y) is precisely the unstable manifold shown in Figure 3. 

We must also fix the function p(e). We expect our results to be sensitive only to 
the e <^ 1 form of p(e) oc poe'^~^. It is thus convenient to choose 



44 



Pie) 



(4.17) 



as a suitably normalized form for p(e). 

Our procedure is now quite simple. Starting at the first branch point, at the root 
of the tree, we fix n, the total number of particles in the tree. We then choose a value 
of e with the distribution given by Eq. (4.17), and compute the values of x and y at 
that branch point from Eqs. (4.3-4). We then repeat the process on each of the two 
sub-branches, and so forth. If the total probability on a branch is zero (as will be the 
case for values of ry > 1 on the weaker branch in model Z), then it will not contribute 
to YliPii be disregarded. If the total mass of a branch is less than one, 

then we stop the subdivision process. We then form the sum either of the probability 
moments of these ends, which correspond to the elementary sub-branches (indexed by 



or else the sum weighted by a simple function of the mass of each elementary sub- 
branch 



We have used the factor n^"^ on the right-hand side of Eq. (4.19) because 0-^(9) = 
aQ{q) = q—1 for a non- fractal cluster, and the elementary sub-branches are presumably 
not fractals. Numerical results for the scaling behavior of these two versions of the 
partition function were, for practical purposes, identical. We can now determine the 
"apparent" dimensions crp(g) by studying the scaling of (logZ(g)) with logn, where 
(• • •) indicates the average over the statistically independent cluster realizations. 

Although we performed simulations for a variety of values of < < 1, we shall 
display our results only for u = 0.6 (results for other values of u are qualitatively 
similar). We used values of n varying from n = 10 to n = 2560, varying by powers of 




(4.18) 




(4.19) 
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two, and we averaged over 10^ realizations at each size. Figure 21 displays (log Z) vs. 
logn for two values of q. The linear scaling is quite evident in these plots. Extracting 
the slope of the linear part, we find (Tp(q') as displayed in Figure 22, with its corre- 
sponding j{a) (defined as the Legendre transform of cr-p{q), not of T(q')-see section 
1 above for a discussion of the relation between these two quantities). Note that the 
f{a) corresponding to the numerical quenched dimensions does not have any negative 
values of / appearing. 

We can implicitly solve Eqs. (4.2) or (4.12) above in order to find (Jj\^{q) or the 
apparent quenched dimensions cr-piQ)- These are displayed with the numerical results 
in Figure 23. For comparison, we also show the one-propagator result from Eq. (2.37). 
It is clear that the annealed dimensions agree with the numerical results only ior q < 2, 
while the apparent quenched dimensions from Eq. (4.12) agree over the larger range 
q < 5. From Eq. (4.14) for nc{q), we already expect that for q = 3 the crossover from 
apparent quenched to annealed dimensions (for model Z with = 0.6) in the numerical 
results will occur only for n ^ 10^. Since these dimensions are essentially identical for 
g < 3, it is not feasible to see this crossover in the numerical results. 

It is possible to understand the fact that for q not too large, the cr-p{q) ~ (^^(g), 
by recasting the expansion in the formal parameter S as an expansion in g — 1. In this 
way, one can show that a-p{q) = cr^(g) + 0{{q — 1)^). This is discussed in detail in 
appendix E. 

5. Conclusions 

We now summarize the main results of this work, which concern the partition 
function Z{q, cr; n) defined in Eq (2.2). 

1. The multifractal dimensions for the branched growth model may be defined ei- 
ther by annealed averaging, for which lim„_»oo(-2(g, cr^(g))) = 0, or by quenched 
averaging, limn^oo{^ogZ{q,aQ{q))) = 0. 

2. It is possible to compute (log 2^) perturbatively by expanding in a formal param- 
eter S, which counts branching vertices in a diagrammatic approach. The result 
is a formal expression 
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oo oo N 

(logZ) = 5^(log2)^ = E E /3^,Mlog^n, (5.1) 

N=l N=l M=0 

where the parameter 5 = 1 in the physical case and the coefficients {Pn,m} can be 
determined from the unstable manifold for branch competition in the underlying 
model. 

3. It is possible to resum the leading M = N and sub- leading M = N — 1 terms at 
all orders N in Eq. (5.1). The result is an expression for (log 2) which combines 
a logarithmic term in n with a constant and a decaying power law: 

(log Z{q, a, n)) = ((^'(g, <T)))To{q, a) logn + ri(g, - n"^) + • • • , (5.2) 
with A(g) > 0. 

4. From Eq. (5.2), we see that in the limit n — > oo, crQ(g) is fixed by requiring 
((•0)) = 0, which is also the criterion that fixes a"_4(g). Thus, in the limit of 
large n, (JQ^q) = cT^(g). This limit is attained only for n Uc = exp(aexp6g). 
For n <C Uc, the apparent multifractal dimensions cr-piq) ^ Cyi(g). Furthermore, 
there is numerical evidence that, unlike the annealed dimensions, these apparent 
dimensions do not exhibit negative values of f{a). 

5. For appreciable values of g — 1, tt-c is enormous, so that the true quenched (i.e., 
annealed) dimensions are virtually unobservable. For g — 1 ~ 1, the true quenched 
dimensions will be observed; but in this case cr-p(g) — > cr^(g) anyway. 

The multifractality of these models is of a quite interesting kind. The average of 
the partition function displays true scaling over the entire range of g. The average of 
its logarithm, on the other hand, contains regions above and below ric with different 
scaling properties. As a function of q, there is a weak function qdn) such that for 
q < qc and q > qc essentially different scaling properties are being explored. Since 
the essence of multifractality is the smooth variation of scaling properties with g, this 
represents, in a subtle sense, a failure of multifractality. 
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It is therefore most accurate to say that these branched growth models are mul- 
tifractal only as an ensemble; typical members of the ensemble exhibit this weak de- 
viation from multifractality. Since the branched growth model seems to be a quite 
adequate theory, both qualitatively and quantitatively, for diffusion-limited aggrega- 
tion, we expect these quite novel properties to hold also for DLA. 
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Appendix A: "Propagator" sums 

It was claimed above in section 2 that 

N-l N 



j>l 3 M=0 



where the sum over j is down a "main hne" with a total of n particles. First we repeat 
the computation of Ai,i. From Eq. (A.l) 



5^^) = Ai,ilogn + Ai,o. (A2) 

>i 



3^ 



Separating out the first term in the sum, we have 



Ai,ilogn + Ai,o = ^ + (Ai,ilog[(l-y)n]+Ai,o), (A3) 



which simplifies to 



^ + (Ai,ilog(l-2/)) = 0, (A4) 

?7' 



or (recalling Eqs. (2.7), (2.8)) 



Ai,i = A = -— — (^.5) 

((log(l-y))) 



The factor of po has cancelled against that coming from the integral, so that Ai,i is 
independent of po; as are all of the {Ajv,m}- The constant Ai^o (elsewhere called oq) is 
not determined by this argument; in fact we expect it to depend on the specific cut-off 
procedure in the theory, and thus be non-universal. 

To compute the remaining {Xn,m}i we simply generalize the above argument. 
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, N-l N 

^1^ :^ — ^) = 2^ ^^'^ " 

j>l j M=0 



Polog " , / X 1 Mr/1 ^ 1\ 

+ ( 2^ >^N,M log [(1 - y)n]) 



M=0 



+ 2^ Aat.m ( ^ J log (log (1 - y))^ 

M=0 L=0 ^ ^ 

(^•6) 

leading immediately to 



M=l L=0 ^ ^ 

The left hand side of Eq. (A. 7) is a sum of terms of the form log^ n, with < 
K < N — 1. Setting each of the coefficients ax equal to zero yields N simultaneous 
linear equations for the N quantities {Ajv,M5-^ ^1}- The solution of these equations 
then gives {Ajv,M,Af > 1}. The non-universal coefficients {Ajv^o} do not appear in 
Eq. (A. 7), and are left un-determined by this procedure. 

As an example, consider Xn,n- These may be obtained from the coefficient of the 
log''^""^ n term in Eq. (A. 7), 

It is also simple to determine the {Ajv,jv-i}- The coefficient of the log^~^ n term in 
Eq. (A.7) is 



Aiv.iv-i (aog(i - y))) + Xn,n ((iog^(i - y))) = 0, (^-9) 



giving 



Ar ((log^(l-y))) l ((log^(l-y))) _ 
2 ((log(l -I/))) 2 ((log(l - y)))2 
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The next coefficient, that for the log^~^ n term from Eq. (A. 7), gives 



^"'"-^ - - 1> 1 6 ((iog(i - y)r - 4 «iog(i - y)r ^ ■ * ' 



Clearly, one may iterate to generate any of the Ajv m that one needs. 



Appendix B: Series in 5 to fourth order 

In section 2 above, we demonstrated the rules for computing (logZ) order by 
order in 5, and computed the terms of 0{5) and 0(5^). The diagrams corresponding 
to these terms are displayed in Figure 24, the results are 



(logZ)i = (Alogn + ao)(W), {B.l) 



and 



(log 2)2 = (Alogn + ao) 



A((V))((/-log(Y3^)))-^((V') 



{B.2) 



Recall that in our more formal notation, A = Ai,i and the non-universal constant 
ao = Ai^o- In the expressions below, we will find it convenient to mix these two 
notations, as well as to use A' = A2,i. 

At third order in S, the diagrams appearing are shown in Figure 25. The attentive 
reader will recall that each diagram will in general involve a number of terms, due to 
the different ways in which propagator terms can be averaged with vertices /_ or g+. 
A lengthy computation yields the following result for (log 2^) 3: 
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(log Z)3 =(A2,2 log^ n + X2,i log n + A2,o) 

• [2A((/_))((V^)) (A((V^log(l - y))) + X'm) + A((/_))((V^2))] 
+ (Ai,ilogn + Ai,o) 

• {A'M(^((V'))((/-log^-^))) + ((/-log(-^)r 
I \2 1-y 1-y 

+ ((/_ log[y(l - y)]))((V'log(l - y)))) (S.3) 

+ AA'((V'f ((/- log[y(l - y)])) + A((^2))((/_ log(l - y))) 

-A((V^))((/_V'log(-^))) + i((V'^))} 
1-y 3 J 

+ (Ai,i logn + Ai,o){((V'))[2Aao((/-))((01og(l - y))) 
+ 2XX2,o{{f-))m + al{{f-))m] + ao((/-))((V'))}, 

where the reader will recall that A2,2 = A/2, the other propagator coefficients being 
given after Eq. (B.2). 

The non-universal terms at this order are collected at the end of the right-hand 
side of Eq. (B.2). Unlike the second order in 5, at 0{5^) there are non-universal terms 
proportional to logn appearing. It can be shown by explicit computation that up to 
0(5^), all such terms come from decaying power laws in n, which in addition vanish 
when ((/_)) — > 0. Thus these non-universal terms survive neither in the limit n — > oo 
nor in the limit q ^ oo. 

We have also performed the 0(5^) computation; we do not show the diagrams. In 
Eq. (B.4) below, we organize the result into four parts: (universal) terms multiplying 
the third-order propagator A3. 3 log' n + A3^2 log n + A3,i logn + A3.0, universal terms 
multiplying the second-order propagator A2,2log n + A2,i logn + A2,o, non-universal 
terms multiplying this propagator, and terms both universal and non-universal multi- 
plying the first order propagator Ai^i logn-|- Ai^O) which we omit. The universal results 
allow the reader to check our re-summation in section 3. The non- universal terms al- 
low the reader to check our assertion that such terms resum to simple decaying power 
laws in n. 
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(log Z)4 =(A3,3 log^ n + A3,2 log^ n + A3,i logn + A3,o) 

- 2A2((/_f ((v^))(A((v^iog(i - y))) + x'm) - x'iif-rm)] 



+ (A2,2 log^ n + A2,i log n + A2,o) 



2X^m [((/_))((V^))((/_ log(l - y) log 

+ mUf- log ^))((V;log(l - y))) 
- 5((/_))((/_ log(l - y)))((V'log(l - y))) 
+ ((/_))((/_ logy))((V'log(l-y)))" 



y 



1 



+ A^A3,2 

+ a2a' 



'((/-log 



1-2/ 



H(/-log^ 



+ 4((/-))((V))((/- log ^)) - 4((/_f ((V^log(l - y))) 



8((/-)) ((V')) ((/- log(l - y))) - 4AA' ^ ((/_ 



+ A^ 



2((/_))((V'))((/_V'log 
((/■ 



y 



i-y" 

Hog{l-y))) + 



'»((/- log 



y 



'2 ' 



- ((V'logii - y))) [{{f-)){{r)) + nnv-n + mv-)) 

-4((/-»((n((/-log(l-y))) 

- AA' (2((/_))((V^))((V'^)) + 2((/_f ((^2)) 

- >^{{{f-9+))m) + ^ ((/-)) ((V'')) + ((/-)) ((V''))) I 

+ (A2,2 log^ n + A2,i log n + A2,o) 

• " - AX\o{{f-rmmog{l - y))) - 2XaUf-r 

- Ax'x,,o{{f-rmf - 2Xao{{f-rm) 

+ (Ai,i logn + Ai,o) •(■■•)• 



(BA) 



53 

Appendix C: Some useful sums 

Our starting point is the standard exponential sum: 



CO 

5^^=exp(x). (ai) 

L=0 



and a modified version thereof 



Now consider the sum appearing in Eq. (3.10), which has the form 

ii=0i2=l 1=0 m=0 \ ^ / V I I y 

where L = li + 12. It is convenient to add and subtract the term Z2 = in Eq. (C.3), 
and set 

S'{x,y)^S{x,y)-So{x,y) {CA) 

with 



00 L L — li 00 

EEEI 

L=Oli=0 1=0 m=0 

and 



r—n;,— n ;— n rr,=n \ / \ ' ' / 



00 00 ^ 

In order to compute S{x,y) from Eq. (C.5) we first alter the order of summation 

L L-h L L-l 

EE -EE. (c.i) 

Zi=o 1=0 1=0 h=0 

leading immediately to 
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CXD L oo 



mill {m + L + 1)1 



x"'{2y)^. (C.8) 



L=0 1^0 m=0 

To perform the remaining sums, we change variables to N = m + L, n = m + I, and 
rewrite the sums as 



oo L oo oo N n 

EEE-EEE 

L=0 1=0 m=0 N=0 n=0 m=0 



We now have 



(C.9) 



oo N n ^ 

si^,y) = E E E TFTW'"""''^'""'"- 

The remaining sums are now elementary. Summing over m, we obtain 



(C.IO) 



and summing over n, 



Recognizing that 



we obtain our final result: 



oo N ^ / \ n 

N=0 n=0 ' / \ y / 



AT 



(%) 



TV 



1 - 



1 I 1+ - 

2 ^ y 



oo ^ 

E " 

iV=0 



'-AN + 1)1 



exp(x) — 1 



a; 



(C.ll) 



(C.12) 



(ai3) 



5(a;,y) = 



exp(2y) - 1 _ 1 A ^ ^\ exp (2y + 



y-x 



(C.14) 
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A similar but simpler algebra yields the sum So{x,y) from Eq. (C.6), 



Soix,y) = ^^^ (C.15) 
y-x 



so that 



S'{x,y)^{ey-1)(^^-^) (C.16) 

V y-x J 

This result illustrates the necessary tricks for performing the types of sums ap- 
pearing in this paper. To give a further example, the sum appearing in Eq. (3.9) is of 
the form 



°° ~ ~ 'K/L-1\ (M + 1 + 1)1. 



Zi=0 l2=0 mi=0 m2=0 1=0 



h J mil im2 + l + 1)1 (M + L + 2)1' 

(C.17) 



with L = li + I2 and M = mi + m2- Once again, 

00 00 ^2 00 L L—l 

«i=0 l2=0 1=0 L=0 1=0 li=0 

SO that 



s"i.y)= E E t±^'-' Ji:':,%HM + L + 2y. ^'y"- (^-1^) 

mi=0m2=0L=0 Z=0 i v ^ ' ' / V ' / 

In addition to M = mi + m2 and L = li + I2, we take N = M + L, and n = M + I. 
We then rearrange the sums 

oocxdcxdL 00 N n n 

E EEE-EEE E- (^-^o) 

mi=0 m2=0 L=0 Z=0 N=0 n=0 mi=0 M=mi 

we find that the sums can now be performed to yield 
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S"{x,y)='-r- ^ 

2\ x-y 

The other sums appearing in section 3 can all be performed using this repertoire 
of tricks. 

Appendix D: Computation of \og{Z) 

Although, as we have seen, the computation of (logZ) is quite intricate, the 
computation of {Z) is by contrast straightforward. This is quite useful, since we have 
seen that the terms in (logZ) that are oc logn as n — > oo are identical to \og{Z). 

Consider the initial branching point in the cluster. At this point, we can write 
the relation 



(C.21) 



/•CO 

{Z{q, (j; n)) = / dtp{t) {f-{en){Z{q, a; y{tn)n)) + f+{en){Z{q, a; (1 - y(en))ri))} , 
Jo 

{D.l) 

where the averages inside the integral on the right-hand side now do not include the 
first average over e. In Ref. 11, we showed how this integral equation can be expanded 
into differential equations of increasing order, whose solution determines the annealed 
dimensions cr^(q'). For general cr these differential equations have solutions of the form 
{Z) = n'^(^'^). Using this Ansatz in Eq. (D.l), we obtain 



POO 

/ dep{e) {f_{enMen)nr + f+{en) [(1 - y{en))nf} . {D.2) 
Jo 



Factoring out n^, substituting /+ = 1 + ^r^, and then transforming to the rj = en 
variable, we see that 



nOO 

/ drj v''-'{f-{v)y''{v) + (1 + ^+(^))[i - y{v)r - 1} 
^0 



{{f.y'^ + {l+g+){l-yr-l))^0, 

(D.S) 
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which imphcitly determines fx{q, a). We can expand Eq. (D.3) so as to obtain order 
by order in S, the formal parameter that counts the number of powers of /_ or g^. 
If /_ — > and g+ — > 0, we see that — > as well. To implement the expansion, 
we write y'^ = exp(//log7/) and (1 — y)'^ = exp(/ilog(l — y)); we can then expand the 
exponentials order by order in fj, in Eq. (D.3), thereby obtaining solutions for /i to any 
desired order m S. To second order in 5, we can write 



+ ^[log(l -y)+g+ log(l -y) + f- logy] + ^ log'(l - y))) + 0{6^) = 0, (DA) 



which gives 



fj, = 



^ _ ((/-logy + ff+log(l-y))) ^ #))((log'(i - y))) 



{D.5) 



((log(l - y))) L ((log(l - y))) 2((log(l - y)F 

or, recalling the definitions of A and A' (Eqs. (A. 5) and (A. 10)), 

= A#))[l + A((/_ logy + g+ log(l - y))) + A'#)) + 0{6% {D.6) 

which agrees with the term oc logn in Eq. (3.43) above. Note that {{ip)) = implies 
that /X = to any order in S. Clearly terms of higher order in S can be obtained from 
Eq. (D.3) with considerably less labor than by direct perturbative computation. 

Appendix E: 5 expansion as an expansion in — 1 

The perturbation expansion derived above can be formally written in the following 
form: 

(log^) = mF,{q, a; n) + U^))F2{q, a; n) + • • • , (E.l) 



because the last vertex (at the far right) of any diagram is always averaged without 
any powers of logy or log(l — y), and can always be either /_, g^, or a product of 
factors of these two. Thus this last vertex can always be regarded as being of the form 
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((V''')). (This property is satisfied by Eqs. (3.43) and (B.1-4).) Recall the definition of 

^(.; a) = /_(.) + .4.) = ^ + - 1. (^.2) 

At = 1, for cr = we have i/j = 0, which implies that both cr^(g) and the apparent 
(7'p{q) are equal to zero. Thus it is clear that 

V'(g,^p(g);^) = c>(g-i). (e.3) 

Since by definition {{i^{q, ^^{q)))) = exactly, Eq. (E.l) together with (E.3) implies 
that a-p{q) = cr^(g) + 0((g - 1)^), with 

mq,ar{qm = 0{{q-lf), {EA) 

and 



m)=oiiq-in 
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Table caption 

1. This table shows the orders of the terms divergent in logn appearing in the S- 
expansion for (log 2^). A solid cross indicates that a term of that order in 6 and 
logn appears; a dashed cross indicates that while such terms appear in some dia- 
grams, they cancel between diagrams. The leading and sub-leading resummations 
discussed in section 3 correspond to sums down the diagonals of this table. 



Figure captions 

1. A typical two dimensional DLA cluster, grown using off-lattice random walkers. 
There are approximately 35, 000 particles in this cluster (courtesy of M. Leibig). 

2. Multifractal scaling functions f{a) for (a) a typical non-stochastic system, and 
(b) a typical stochastic system, with annealed averaging. Note the appearance of 
negative values of / in the latter case. These functions are obtained by Legendre 
transformation of T(g) or a(q) oc T(g). 

3. The unstable manifold for the dynamics of branch competition in "model Z" . The 
two branches move along the indicated diagonal line until either x = or x = 1. 
This corresponds to the complete screening of one or the other of the two sibling 
branches. From this point only the parameter y reflecting the relative masses of 
the two branches changes. 

4. The unstable manifold for the dynamics of branch competition for DLA in two 
dimensions, as calculated by a renormalization method in Ref. 12. Comparing 
with model Z, we see that screening is a more gradual process in DLA. 

5. Two sibling elementary sub-branches. The last node in the branching cluster is 
indexed by J, has a total of nj descendants, and has a total growth probability 
of pj. The stochastic variable controlling the branching at this node is ej. The 
stronger descendant sub-branch has a number of particles ng = (1 — y{^j'nj))nj, 
and a total growth probability = (1 — x{ejnj))pj 

6. The sub- branch found at the bottom of the main branch, defined by taking at 
every node the stronger of the two paths down the branch. The index j indexes 
the nodes along this main branch. 

7. The terms of first order in 6 containing one factor of /_ arise from the elementary 
sub-branches at the end of the "first-order" sidebranches off of the main branch. 
The terms of first order in S containing one factor of g+ arise from corrections 
to the partition function contribution coming from the elementary sub-branch at 
the end of the main branch. 

8. Diagrams at 0{5) in perturbation theory. The solid circle indicates an /_ vertex, 
and the open circle indicates a g+ vertex. A solid and an open circle connected 



by a short vertical line indicate a vertex of ip = + g+. A horizontal segment 
indicates a branch over which n~'^ is summed, possibly with a logarithmic factor 
in n. 

9. A diagram at 0{5'^) that is oc g'^. This diagram comes from {Z — 1). The two 
horizontal bars indicate the two propagator sums. Note that including a factor of 
Qj^ does not lead one onto a weaker side-branch. 

10. A diagram at that is oc This diagram comes from (Z — 1). The two 
horizontal bars indicate the two propagator sums. The vertical bar indicates that 
the rightmost /_ vertex is off of the main branch. Taking a factor of /_ in a 
diagram arising from {Z — 1) always moves one onto a weaker sidebranch. 

11. Mixed terms (x ■ g+ arising from {Z — 1). These diagrams are 0(5^). 

12. All terms of 0(5^) arising from — (1/2)((Z — 1)^). Note the symmetry factors. 

13. All diagrams containing only one propagator in the full expression for (log 2^). 
Note that the joined vertices indicate factors of /_ or g^ averaged together, 
e.g. ((/_(7+)) for the last diagram. The sum of these diagrams gives the "one- 
propagator" approximation {\ogZ)ip. 

14. A typical contribution to {Z — 1). At leading order, the leading logarithm in each 
propagator is taken. Each vertical segment indicates a move to a side-branch. 

15. A typical contribution to {{Z — 1)^), arising from a product of a diagram with li 
factors of g+ upstream from the first factor of /_ , with mi factors of ip off of the 
main branch, with a diagram with I2 factors of g+ upstream from the first factor 
of /_ , with 777-2 factors of i/j off of the main branch. For simplicity, we have drawn 
all factors of t/j along the same branch; terms involving more than one factor of 
/_ will properly have further vertical segments. (Note that wc use "upstream" to 
mean closer to the root, at the left, not closer to the elementary sub- branches.) 

16. One term arising from the product indicated in Figure 15. Defining L = li + 12, 
we have L — I > li factors of g+ upstream of the first factor of /_, followed by 
I < I2 factors of g^^ and then the second factor of /_. The diagram includes mi 
factors of ij) coming off of the first /_ vertex, and m2 factors of ij) coming off of 
the second /_ vertex. 

17. This diagram indicates a term in the product indicated in Figure 15 that is oc 



{{f-))- Now all L = I1 + I2 factors of g+ are to the left of the vertex containing the 
two factors of /_ averaged together. Note that at the sub-leading order, which 
we are considering here, the two separate side-branches emanating from these two 
factors of /_ (one containing mi factors of i/j, and the other containing 777.2 factors 
of V') can be treated as independent. 

18. This diagram indicates a term oc ((/^)) arising from {{Z — 1)^)). At sub-leading 
order, the factor of 2 — 1 whose first /_ vertex appears separately can be averaged 
independently of the other two factors of Z — 1. 

19. At sub-leading order, terms oc ((/_logy)) can be obtained by averaging a single 
factor of Z — 1; a typical diagram is shown. The arrows indicate the /_ vertices, 
each of which is a potential source of a factor of ((/_ logy)) except the last, which 
has no propagator to the right to provide a factor of log y. 

20. At sub- leading order, terms oc ((/_ log(l — y))) require at least two factors of Z — 1. 
This diagram arises from a product of the form shown. The arrow in the diagram 
corresponding to a single factor of Z — 1 indicates the vertex (with a factor of one 
rather than a factor of /_ or g^) at which the term log(l — y), arising from the 
propagator, is taken. To the right of this arrow, it is not necessary to keep track 
of whether vertices are /_ or g^, provided at least one such vertex appears. In 
the combination of the two factors of Z — 1, this log(l — y) is averaged with a 
factor of /_ contributed by the other partition function. 

21. Numerical results for (logj^^p^) vs. logn for a Monte Carlo realization of model 
Z random branched growth, for q = 2 and q = 5. We used a parameter 1/ = 0.6 
to specify model Z. The linear slope indicates apparent multifractality, the value 
of the slope giving the apparent dimensions a-piq). 

22. (a) The apparent dimensions <J-p{q) vs. q for Q < q < 10, and (b) the Legendre 
transform of this function, /(ck), computed numerically for model Z with 1/ = 0.6. 
No negative values of / appear. In addition, since we compute cr>i.(g) only for 

> 0, we show only the left side of the multifractal spectrum. Note that model 
Z has the pathological feature that limq_»o (T{q) ^ —1, due to the fact that much 
of the cluster surface has growth probability zero. 

23. cr vs. q for model Z with v = 0.6 computed in various manners, (a) solid line: 



numerical results, as in Figure 22(a); (b) dashed line: annealed dimensions cr^(g); 
these are also the apparent dimensions a-p{q) at 0{5); (c) dot-dashed line: ap- 
parent dimensions, computed to 0(6'^), and (d) dotted line: the one-propagator 
result for a{q), from Eq. (2.37). We note that the computation for apparent di- 
mensions gives excellent results ior q < 5, but fails at higher values of q. The 
one-propagator result gives qualitatively correct results over the entire range of q. 

24. Diagrams at 0{6) and 0(5^). These correspond to the terms appearing in the 
text in Eqs. (B.1-2). They also summarize Figs. (8-12). 

25. Diagrams at 0{5^). These correspond to the terms appearing in Eq. (B.3). 
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